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Abstract 

The unpolarized response functions of the quasielastic ie O(e, e'p) 15 N reac- 
tion are calculated for three different types of relativistic bound state wave 
functions. The wave functions are obtained from relativistic Hartree, rela- 
tivistic Hartree-Fock and density dependent relativistic Hartree calculations 
that reproduce the experimental rms charge radius of 16 0. The sensitivity 
of the unpolarized response functions to the single particle structure of the 
different models is investigated in the relativistic plane wave impulse approx- 
imation. Redistributions of the momentum dependence in the longitudinal 
and transverse response function can be related to the binding energy of the 
single particle states. The interference responses Rlt and Rtt reveal a strong 
sensitivity to the small component of the relativistic bound state wave func- 
tion. 

PACS number(s): 21.30.Fe, 21.60.Jz, 24.10.Jv, 25.30.Fj 
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I. INTRODUCTION 



The exclusive (e, e'p) proton knockout reaction is a powerful tool for the investigation of 
the single particle structure of complex nuclei [jl]-§J. Quasiparticle properties, such as oc- 
cupation probabilities, spectroscopic factors, binding energies and momentum distributions 
can be determined and compared to the results of theoretical models, which have to include 
both, information about the electromagnetic reaction mechanism and the wave function of 
the bound nucleon. 

The nuclear structure contributions can be seen e.g. in the reduced cross section. In the 
nonrelativistic plane wave impulse approximation, the reduced cross section is proportional 
to the momentum distribution of a single bound proton inside the nucleus. This propor- 
tionality is characterized by the spectroscopic factor, which corresponds to the probability 
that a single particle state is occupied by a nucleon. Therefore, spectroscopic factors can be 
extracted comparing experimental scattering data with theoretical predictions for the mo- 
mentum distribution of the nucleus. In a mean field calculation the spectroscopic factor is 
equal to one for occupied or zero for unoccupied states, respectively. Different spectroscopic 
factors indicate the deviation from the shell model picture and therefore the importance of 
nucleon-nucleon correlations. 

The theoretical background of the scattering formalism is provided by quantum electrody- 
namics, a complete relativistic framework, which describes the electromagnetic interaction 
with highest accuracy. Consequently, all contributions to the scattering amplitude espe- 
cially the hadronic current operator are relativistic expressions and the matrix elements of 
the current should be calculated between states obtained from a relativistic treatment of the 
many-body problem. Due to the higher complexity of the relativistic problem the solution 
of a nonrelativistic reduction is more economic and it is able to reproduce the experimental 
data in a wide range of missing energies and momenta. 

Relativistic calculations of exclusive (e, e'p) scattering reactions, including electron distortion 
and final state interactions, have been performed @-f§. The results provided spectroscopic 
factors e.g. for the 3si/ 2 and 2di/ 2 shells in 208 Pb of S a ~ 0.7, consistent with earlier theo- 
retical predictions. Though these spectroscopic factors are extracted from the low p m data 
(Pm < 300Me\^), relativistic calculations can give simultaneously a good reproduction of 
the high p m data 0, where the main effects arise from the improved relativistic treatment 
of the electron distortion and the final state interaction. 

The inclusion of electron distortion and final state interactions can give a good descrip- 
tion of the experimental data, but for a deeper understanding of the relativistic reaction 
mechanisms the relativistic plane wave impuls approximation (RPWIA) seems to be an ap- 
propriate calculation scheme. Recently, the RPWIA was chosen to study the role of the 
negative energy components of the bound nucleon wave functions || . Analyzing the factor- 
ization of the scattering cross section, a feature of the nonrelativistic limit, the importance 
of relativistic effects for different choices of current operators, kinematics and restorations 
of current conservation has been investigated thoroughly. All calculations were made using 
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the bound state wave function of Horowitz and Serot f9HTT]|. Therefore it is interesting to 
ask, what influence on the results can be observed, if bound state wave functions from other 
relativistic nuclear structure calculations are considered. This question shall be studied in 
this paper. 

For the investigation bound state wave functions from three different approaches were chosen, 
namely, the Relativistic Hartree (RH) approach of Horowitz and Serot |9Hill , the Relativis- 



tic Hartree-Fock calculations (RHF) of Bouyssy et al. |12] and a Density Dependent Hartree 
(RDDH) approach including rearrangement terms of Fuchs and Lenske |I3| . Each of them 
was used to calculate the response functions of the 16 0(e, e'p) l6 N electron scattering reac- 
tion. 

All models used are based on a microscopic understanding of the nucleus using neutrons 
and protons as effective degrees of freedom and the exchange of a-, u>-, p- and, in the RHF 
approximation, 7r-mesons to mediate the nuclear force. They provide a consistent mathemat- 
ical description starting from a covariant Lagrangian. In these models the most important 
contributions to the nucleon-nucleon potential arise from an attractive a-meson exchange, 
which is understood as a parametrization of the 2tc exchange diagrams, and a repulsive 
o;-meson. Calculations in finite nuclei within the mean field approximation are character- 
ized by two large potentials of scalar and vector type (S-V) cancelling each other to a large 
amount. As a result the spin orbit splitting emerges automatically. 

The RH and the RHF approaches use phenomenological one boson potentials. In both 
models the model parameters are determined in essentially the same way. The er-N and 
cj-N coupling constants and the cx-meson mass were fixed in both cases to reproduce the 
saturation point in nuclear matter and the charge rms radius of 16 O. 

The third model is the RDDH approximation, which is a first step to the relativistic 
description of finite nuclei using realistic forces. In the RDDH approximation the rela- 
tivistic Brueckner-Hartree-Fock (RBHF) potential of a realistic interaction (G matrix) is 
parametrized in nuclear matter in terms of density dependent coupling constants and the 
parametrized interaction is then applied to finite systems. This calculation scheme can 
be regarded as a reliable approximation for the self consistent solution of the relativistic 
Brueckner-Hartree-Fock equations in finite nuclei RDDH calculations were performed 
by Brockmann and Toki [T5]. It was extended to the Hartree-Fock approximation by Fritz 
and Muther [IB] and Boersma and Malfliet [|T7|. In the work presented here, we use bound 



state wave functions from the vector density dependent RDDH approach of Ref. [|13|]. This 



many body approximation accounts for additional rearrangement terms arising from the 
field theoretical argument that the density dependent coupling constants have to be consid- 
ered as functionals of the baryon fields. Due to the rearrangement terms this model goes 
beyond an effective Brueckner Hartree-Fock approximation. The model was chosen from 
other RDDH calculations due to its good results for the single particle as well as the bulk 
properties of 16 O. 

The results for the three nuclear structure calculations yield essentially identical results for 
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the charge radii of 16 0, but have a different single particle structure. The single particle 
structure can be tested directly with (e, e'p) knockout reactions. Using this reaction we 
want to investigate the sensitivity of the unpolarized response functions in the RPWIA to 
the model functions and to relativistic effects included in the bound state wave functions of 
the three approaches. 



II. ELECTRON SCATTERING OBSERVABLES 

As stated in section 0, the theoretical description has to account for the scattering for- 
malism as well as for the description of finite nucleon systems. In the following section the 
theory will be summarized, which is needed for the electron scattering observables calculated 
in the present work. The theoretical foundations for exclusive electron scattering reactions 
were developed in several publications |18 |- p5| . The present paper follows the conventions 



of |25] except to some changes in the notation. The nuclear information is included in the 



matrix elements of the hadronic current 

J^q) = (p,s x ^ fPB \J^q)\^p A ) . (1) 

The initial state ipip A describes the A particle target nucleus with the four momentum 
P% = (Ma, 0), whereas the final state consists of the A — 1 particle residual nucleus ipfp B 
(Pg = (Eb,Pb)) an d the ejected proton (p% = (E x ,p x )). The present investigation will be 
performed in the RPWIA, i.e., the outgoing proton is represented by a plane wave u(p x , s x ) 
and the current is given by the current of an individual constituent, treated as a free particle. 
The full complexity of the initial and outgoing states for the bound nucleon system is 
therefore replaced by a single shell model wave function ip a (p). In this approximation one 
can write for the above current matrix elements 

J»(q)=u(p x ,s x )J»{q)iP a {p)- (2) 

The calculation is done completely in momentum space. The bound state wave function 
ip a {p) for the initial state is the information we obtain from the shell model calculations 
and p = p m corresponds to the missing momentum. The results of the approximation, 
neglecting proton and electron distortion, should not be compared to the experimental data. 
This approximation, however, should allow us to study the sensitivity of the results on the 
model for the bound state wave function ip a (p). 

In the present investigation the main attention will be paid to the response functions, which 
are directly related to the hadronic tensor. If polarization is not taken into account for the 
incoming electron beam, as well as for the target and the outgoing particles, we can define 
the hadronic tensor including an average over the initial states and a sum over the final 
states 

w" = £ $ -a-") •/*"(?) J v (q) • (3) 
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In this definition an integration over the momentum conserving delta function has already 
been performed, which provides the relation q = p x — p. The response functions can be 
constructed using the hadronic tensor according to 



(4) 



Rl = 


/. de p* 

J line 




Rt 


/ de Px 

J line 


(W n + W 22 ) . 


1 > : COS = 


/ de Px 

J line 


(-W 01 - w w 


Rtt cos 20 = 


/ de Px 

J line 





including an integration over a linewidth in the missing energy spectrum. These response 
functions can be used to express the unpolarized scattering cross section as a contraction of 
hadronic responses and the appropriate electron contributions, which are defined as in PJ2~T 



da m \p x \ 



dek>d£lk>d£lp x (2ir) 



a M (V L R L + V T R T + V LT R LT cos + V tt R T t cos 20) . (5) 



Since recoil effects of the residual nucleus will not be considered in the present work, an ap- 
propriate factor has been neglected in the above formula. Using the definitions of this section 
the observables depend, except to the bound state wave function, on standard expressions 
like the current operator or the Dirac spinor. 



III. RELATIVISTIC BOUND STATE WAVE FUNCTIONS 

It was shown in the last section that the information for the calculation of the current 
matrix elements, which is required from relativistic nuclear structure calculations, is the 
momentum space wave function ^f a (p) of a particular shell model state a. The solutions 
of a relativistic Hartree or Hartree-Fock calculation in finite nuclei are usually given in 
coordinate space. The corresponding momentum space wave function can be obtained by a 
Fourier transformation according to 

**(P) = (^372 l d * r e-^Vair) . (6) 

More details about the bound state wave functions are given in Appendix A. In the context of 
(e, e'p) scattering calculations a nuclear structure model should provide a good reproduction 
of the charge distribution, which can be observed directly in electron scattering experiments 
at lower energies. The three nuclear structure models chosen for the investigation here 
satisfy this boundary condition. We want to study the influence of the different single 
particle structure of these models on the unpolarized response functions of Eq. (f|). 
The models we use for the study are the relativistic Hartree approximation (RH) of Ref. [ID 



the relativistic Hartree-Fock (RHF) approximation of Ref. ]TI| and the density dependent 
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Hartree (RDDH) approximation of Ref. [13 . 



The RH and the RHF models follow the same way to fix their free parameters. The cr-N 
and o>-N coupling constants and the cx-meson mass are adjusted to reproduce both, the 
saturation point of nuclear matter and the charge rms radius of 16 0. 

In the Hartree- Fock approximation the Dirac structure of the nucleon self energy S(r) is 
modified due to the nonlocal structure of the two particle interaction. The self energy £(r), 
which incorporates the influence of all other nucleons in the nucleus on a single particle, is 
included in the free Dirac equation according to 

(i0 M y*-m-£(r))tf a (r) = O. (7) 

In a nuclear system, which is characterised by rotational invariance, parity conservation and 
time reversal invariance, the general structure of self energy is given by 

S(r) = S s (r)-7oS (r) + 7-S,,(r). (8) 

In the Hartree approximation the spatial part of the vector self energy S„ vanishes. The 
remaining structure is normally considered as a potential of scalar and vector (S-V) type. 
In the Hartree-Fock approximation the additional X„ term arises and gives a structural 
modification compared to the mean field approximation. 

In the Hartree-Fock calculations the potentials are not given in the separated form of Eq. 
(§). In principle, the Dirac structure can be projected out if the Fock contributions are rede- 
fined in terms of local single particle potentials |T2j and incorporated in the Dirac equation. 
From the nuclear matter calculations we know, that in the Hartree-Fock approximation, 
every meson contributes to E s , E and £„, and that the H v term is small. This should 
remain true also in finite systems. 

The models for the Hartree and the Hartree-Fock approximation are described by the fol- 
lowing parameter sets: The a-meson mass, especially sensitive to the rms radius, is chosen 



as m a = 520 MeV for the RH model |IJ and m a = 440 MeV for the RHF model 0. The 
coupling constants g a and g w were fixed to reproduce the saturation point in nuclear matter. 
In the RH as well as the RHF model the masses of the u- and p-meson are = 783 MeV 
and m p = 770 MeV in accordance with the experimental values. In the RH model [IU| the 
p-meson coupling constant g p /4:iT = 5.19 is fixed to yield a bulk symmetry energy of 35 MeV. 



For the Hartree-Fock calculation we have chosen the parameter set (e) of Ref. ||12|| , which 
provides the best results for the bulk and single particle properties of 16 0. In this set the 
7r-N and p-N coupling constants of the Hartree-Fock calculations were fixed to the physical 
values /^/47r = 0.08 and g p /4:7i = 0.55, where a ratio of f p /g p = 3.7 for the tensor coupling 
of the p-meson was chosen. The 7r-meson mass in this model is m v = 138 MeV. 
For the effective relativistic Brueckner-Hartree-Fock calculation in finite nuclei we have cho- 
sen the vector density dependent (VDD) Hartree approach of Ref. [13]. Here, for the uo- 



meson and p-meson mass as well as for the p-meson coupling constant, the values of fTC 



were adopted. For the a-meson m a = 550 MeV was used. The coupling constants g a and 
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g w were taken from a parametrization of the nuclear matter DBHF self energies [fj6| with 
second order polynomials [^7J . 

The single particle and bulk properties of the three models, which provide the bound state 
wave functions needed in the electron scattering calculations, are summarized in Tab. |. The 
calculated energy per nucleon has been corrected in the RH and RHF models to account 
for the effects of the spurious center of mass motion. A contribution of 0.67 MeV has been 
added to the binding energy. In the RDDH model the results of Ref. fl3|] include already 
the center of mass motion. Since the calculations for the RH and the RHF approach have 
been performed with the computer code developed by Fritz and Miither [ 16| , the results of 
Tab. | differ slightly from the original values given in the Refs. f[0,12| due to numerical un- 
certainties. In the calculation scheme of Ref. |16j the Hartree or Hartree-Fock Hamiltonian 
is diagonalized in a box basis. 



IV. BOUND STATES AND RESPONSE FUNCTIONS 

The reaction considered in the present investigation is the proton knockout reaction 
16 0(e, e'p) 15 N, where the residual nucleus is left in its ground state, i.e. the proton is knocked 
out of the pi/2 state of 16 0. 

To calculate the response functions of the (e, e'p) reaction the bound state wave functions 
of the three models discussed above were used. The computations have been performed 
completely in momentum space. A new computer code has been developed. Note, that for 
the definition of our interference response functions in Eq. (|4]) we use a different convention. 
For the calculation of the knockout reaction we have chosen perpendicular kinematics, which 
provides us with the interference terms of the relativistic response functions Rlt and Rtt- 



For the current operator we take the CC1 description P8[ , using the common dipol form 
factors [Q. It was shown that this operator gives a simultaneous good description for high 
as well as low missing momenta in 208 Pb. Furthermore, this operator enhances the 
contributions from the negative energy projections of the bound state wave functions as 
stated in Ref. ||. For the four momentum transfer g M = — k ffJ- = (uj, q) we take kinematics 
I of Ref. § with 

\q\ = 500 MeV/c uj = 131.56 MeV , (9) 

where the value of uj corrsponds to the quasielastic peak value ujqe- Note, that for the 
definition of our interference response functions in Eq. we use a different convention. 
The response functions for the three models presented in the last section are displayed in 
Fig. [I]. On the upper left panel we see the longitudinal response function R^. In all figures 
the results obtained from RH bound state wave functions are displayed with solid lines, 
whereas the RHF and the RDDH results are displayed with dash-dotted and dashed lines, 
respectively. It can be observed that the maximal strength of Rl is reduced by the RHF 
and the RDDH approach in the order of 10% relative to the RH approach. 
The integrated strength of the response functions 
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h = I dp m p 2 m Ri(p m ) , 



(10) 



(i = L, T, LT, TT) shows if either a reduction corresponds to a redistribution of the strength 
in momentum space or there occurs a true reduction of the strength. The results for the 
integrations are displayed in Tab. P]. For the longitudinal response function one finds values 
of II = 2.81 for the RDDH and II = 2.83 in the case of the RH and the RHF approach, 
where the factor 100 has been multiplied to the original values. Therefore, the reduction of 
the maxima can be understood essentially as a redistribution of the total strength over the 
momentum scale, leading to higher values in the RHF and RDDH response functions for 
large p m . 

For the transversal response Rt displayed in the upper right panel of Fig. [l| a similar be- 
haviour as for Rl is found considering the RH and the RHF bound state wave functions. 
However, the response function Rt obtained from the RDDH bound state wave functions 
shows a stronger reduction relative to the other transversal responses. This becomes more 
obvious, if we consider the integrated strength, which takes It = 3.30 for the RH and 
It = 3.31 for the RHF response, whereas for the RDDH approach we find It = 3.12. 
In the interference contributions this integral reduction can be observed as well. In addition 
a small reduction is induced by the RH wave functions compared to the response Rlt of 
the RHF wave functions, which is reflected by a smaller value of Ilt- This reduction is even 
enhanced in the Rtt response. Note, that the maximum values arising from the RH wave 
functions are larger. That means we have an additional enhancement of the large momen- 
tum components in the RHF transverse responses arising from the increased total strength, 
which is not so obvious in Fig. [I]. 

In order to summarize these features, one may say that we observe two effects 

• a redistribution of the strength of the response functions from lower momenta to 
high missing momenta if we compare results derived from RH, RHF and RDDH, 
respectively. 

• a reduction of the integrated strength if we compare the results of RH and RHF on 
one side with these derived from RDDH on the other side. This reduction is negligible 
for the longitudinal response Rl, but significant for the response functions Rlt and 



In order to understand the redistribution of the strength in the example of the longitudinal 
response Rl-, we consider the single particle densities of the three models in the momen- 
tum and the coordinate space. The single particle density is defined in momentum space 
according to 




(11) 



2jg + 1 
47T 



{dip) + faiP)) , 
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where the upper and lower components of the relativistic wave functions, g a (p) and f a (p), 
are defined in Appendix A. The density in coordinate space is defined in analogy to Eq. ( [TT]) 
with the Fourier transformed expressions. The densities are normalized to the number of 
particles in a particular shell a. 

On the left panel of Fig. |2| the density in coordinate space multiplied with the square of the 
radius, r 2 n pl / 2 (r), is shown. For this expression the area below the curves is identical for 
the three different models. The density n p i/2(r) itself shows an analog but inverse behaviour 
as the longitudinal response function R^. The function r 2 n p i/ 2 (r) indicates that the single 
particle densities of the RHF and the RDDH approaches are shifted to smaller radii. This 
is the reason of the higher momentum components in the wave functions. Though the shift 
in the radial wave functions is small, we have seen that this leads to a 10% reduction of the 
maxima of Rl. This redistribution of the momentum distribution and density distribution 
in configuration space is also connected to the single particle energy of the pi/ 2 state (see 
table |). This state is most weakly bound in the RH approximation. Therefore it is less 
localized in configuration space and the momentum distribution of this state is shifted to 
smaller momenta in this approach as compared to the other calculations. 
In the nonrelativistic PWIA the momentum density of a single particle state is proportional 
to the reduced cross section. On the right panel of Fig. |2] the momentum density n p i/2(p m ) 
is shown in a logarithmic scale, which emphasizes the high momentum contributions. Corre- 
sponding to the longitudinal response function the maxima of n p i/2(p m ) are reduced starting 
from n™£% = 83.3 for the RH wave functions to n™£% = 74.8 respectively n™-/% = 73.0 for 
the RHF and the RDDH wave functions (all values in [(GeV/c)~3\). The numerical values 
show that the behaviour of the momentum density explains the structure of the response 
function Rl. In Fig. § we can see in addition that the upper component of the bound state 
wave function g p i/2(p m ) already reflects the typical momentum distribution of Rl- 
In the RH and the RHF responses beside the effects arising from the different single parti- 
cle densities no significant differences can be observed. The S„ term, included in the RHF 
potential, can therefore be expected to give no substantial modification of the response func- 
tions. 

The reduction of the integrated RDDH responses Rt, Rlt and Rtt has another origin. This 
reduction is closely related to the lower component of the bound state wave function f a . In 
relativistic nuclear structure calculations this small component is enhanced as compared to 
a small component which just arises from a boost of a nucleon with mass m. In order to 
explore the sensitivity we compare results of the relativistic calculation, with those in which 
this enhancement of the small component is suppressed. This nonrelativistic limit \l/™ r (p) 
of the bound state wave function is defined as 

9*(P) = N nr g Q (p) f?(p) = N^-P— g a {p) , (12) 

h p + TTl 

where N nr is a normalisation constant and E p = y/p 2 + m 2 . The small component /™ r (p) is 
just a result of boosting the single particle spinor and therefore of pure kinematical origin. 
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This definition (|T|) means that the momentum dependence of the lower component fa r {p) 
is given, except to a relativistic factor, by the upper component g a (p) of the relativistic 
wave function. On the right panel of Fig. |3] the upper component of this nonrelativistic 
limit is shown. The relativistic upper component 5^1/2 (p) is slightly smaller than than the 
nonrelativistic gp[/ 2 (p) due to the normalisation constants N nr . The situation is completely 
different for the lower component f p i/ 2 (p), which is displayed on the left panel of Fig. f|. We 
find an enhancement of approximately 60% compared to the nonrelativistic wave function 
fp{/ 2 (p) shown on the right panel of Fig. f|. This enhancement of the lower component for 
relativistic nucleons in the medium is consistent with the results found in nuclear matter, 
where the enhancement of the lower component of the Dirac spinor is characterized by an 
effective mass of approximately m* = m + Us = 600 MeV compared to m = 938.9 MeV for 
free nucleons. This additional medium dependence is typical for relativistic calculations and 
it is not included in a nonrelativistic reduction. 

For the relativistic lower component fpi/^ip), the bound state wave function of the RDDH 
approach differs significantly from the lower components of the two other approaches. This 
indicates that the scalar potential Us of the RDDH calculations is less attractive. We have 
seen that in the interference response functions Rlt and Rtt, shown in Fig. [I], these differ- 
ences modify the results notably, whereas the transversal response function shows small 
changes, since it is dominated by the upper component of the bound state wave function. 
The sensitivity of the interference responses on relativistic effects in the bound state wave 
functions can also be demonstrated if the response functions are calculated in the nonrel- 
ativistic limit. This is shown in Fig. A strong reduction especially in the reponses Rlt 
and Rtt can be observed. The strong dependence on the lower component has already 
been shown in Ref. |§. In addition, we can see here that the model dependent information, 
included in the lower component, disappears. The differences of the RDDH responses com- 
pared to the other approaches are strongly reduced. The same effect as discussed above can 
be observed in the nonrelativistic limit of the RDDH lower component fp{/ 2 {p) displayed in 
Fig. [|. In this case the lower components now simply reflect the behaviour of the upper 
components. 

In order to learn more about the typical structure of bound state wave functions, the cal- 
culations have also been performed with the density dependent models of Fritz and Miither 
r5fl . This calculation scheme provides a good description of the binding energies, whereas 



the charge radii are too small. The wave functions are therefore inadequate for a study in 
the context of (e, e'p) reactions and the results are not presented in detail. Nevertheless, 
two additional informations can be obtained from these calculations. The same structural 
difference of the longitudinal response Rl has been found between the Hartree and the 
Hartree-Fock approximation. A reduction of the maxima of approximately 10% and larger 
high momentum components. This can also be traced back to a different single particle 
structure. Furthermore, no reduction in the lower component of the wave function can been 
observed, i.e., the reduction in the RDDH approach of Ref. [O is due to the rearrangement 
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terms and not part of the density dependent approach. 

V. SUMMARY AND CONCLUSIONS 

Bound state wave functions obtained from a relativistic Hartree, Hartree-Fock and a 
density dependent Hartree approach were compared and the influence on the unpolarized 
response functions of the 16 0(e, e'p) 15 N reaction was studied. The nuclear structure models 
chosen reproduce the experimental charge radius of 16 0. It was therefore possible to study 
the sensitivity of the response functions to the different single particle structure provided by 
the models. The aim of this investigation is to give an estimate of possible uncertainities in 
the response functions and to show the influence of the additional information provided by 
a relativistic treatment of the many body problem. 

The calculations were performed in the relativistic plane wave impulse approximation, where 
modifications of the results can clearly be assigned to the bound state wave functions. 
Comparing the three models the maxima of the Rl and Rt responses differ by 10% — 
15%. This reduction of the RHF and the RDDH response functions is accompanied by an 
enhancement at larger missing momenta, so that the integrated strength remains unchanged. 
This behaviour can be connected with the density distribution of the single particle states 
reflecting the single particle binding energies. The vector self energy S„, which appears in 
the Hartree-Fock approximation, has no significant influence on the structure of the response 
functions. This result was expected, since in nuclear matter calculations its contributions 
are found to be small. 

The interference responses Rlt and Rtt are known to be sensitive to the lower component 
of the relativistic bound state wave function. The RDDH lower component is significantly 
reduced compared to the wave functions of the other models. This reduction is clearly 
measured by the interference responses and influences also the integrated strengths of the 
response functions. The maxima of Rlt and Rtt even show a variation of around 20%. In 
the nonrelativistic limit this model dependent information is lost and only small differences 
for the interference responses can be observed. Especially, the interference responses are 
therefore sensitive to a complete relativistic description of the nuclear structure problem 
and can be used to test different relativistic approximation schemes of nuclear many body 
systems. 
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APPENDIX A 

This section gives a survey of the explicit definitions for the relativistic bound state wave 
functions used in the present work. The solutions of the relativistic Hartree or Hartree-Fock 
equations are calculated in coordinate space and have the general structure 

a = {n(ls)jmjTm T } is an abbreviation for the various quantum numbers, where n denotes 
the radial quantum number, / the angular momentum, s the spin, j the coupled angular 
momentum, r the isospin, rrij and m T the orientation of the coupled angular momentum 
and the isopin. The angular momentum of the lower component I' is given by the relation 

l'=\ l + l f ° r l = j-V2l (M) 
1 for l=j + l/2 J U4j 

The normalization of the states is 

Jd 3 r^Ur)^ a (r)= /Vdr (g 2 a (r) + f a (r)) = 1 . (15) 
The Fourier transformation of the relativistic bound nucleon wave functions is given by 

*a(p) = 7^72 j*r e -*""* a (r) (16) 

which leads to the following explicit form of the bound state wave function in momentum 
space 



where the real amplitudes g a {p) and f a {p) are defined according to 



^^K^BmwF'' (17) 
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- J Q r 2 dr ji(pr)gZ-(r) 

fni'M = ]fl fr'dr Mpr) f^(r) (18) 

and for the normalisation of the states it follows that 

J d 3 p^ a (p)^a(p) = J o °°p 2 dp (g 2 a (p) + flip)) = 1 . (19) 

These momentum wave functions served as the input for the relativistic quasielastic electron 
scattering calculations. All calculations are performed in momentum space. Therefore, the 
momentum space wave functions were multiplied directly with the current operator and the 
Dirac spinor to obtain the matrix elements of the hadronic current. 
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TABLES 



RH 


16Q 

RHF 


RDDH 


Exp. 


ls 1/2 [MeV] -37.2 
lp 3/2 [MeV] -16.7 
lp 1/2 [MeV] -8.8 
E/A [MeV] -5.56 
i? cft [/m] 2.75 


-38.6 
-18.0 
-10.6 
-6.10 
2.74 




-36.0 
-17.1 
-12.9 
-7.82 
2.75 


-40±8 
-18.4 
-12.1 
-7.98 
2.73 


TABLE I. The proton single particle energies [MeV] for the lsi/2, IP3/2 an d IP1/2 shells, the 
energy per nucleon (E/A) [MeV] and the rms charge radius Rch [f m ] of 16 for the relativistic 
Hartree (RH), the relativistic Hartree-Fock (RHF) and the relativistic density dependent Hartree 
(RDDH) approach. The calculated energy has been corrected to account for the effects of the 
spurious center of mass motion. 


RH 


16Q 


RHF 




RDDH 


I L 2.83 
It 3.30 
I LT -2.43 
Itt 0.28 




2.83 
3.31 
-2.48 
0.29 




2.81 
3.12 
-2.13 
0.24 


TABLE II. The integrated strength Ii of the various response functions in dimensionless vari- 
ables. The values calculated according to Eq. (10) were multiplied by the factor 100. 
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0.10 
0.06 
0.02 

-0.02 
-0.04 
-0.06 
-0.08 



1 \ ~ 


l 1 l 1 l 

- /c\ RH : 

- if \ RHF 

- 1/ V RDDH - 

; / r t \ ; 


- \ / 

i.i.i 


- — i — i — ' — i — ' — 

Rtt 

- ' X 



0.10 
0.06 
0.02 

0.006 
0.004 
0.002 



100 200 300 100 200 300 
P m [MeV/c] 

FIG. 1. Response functions Rl, Rt, Rlt and Rtt of Eq. (Q) (in [/m 3 ]) as a function of the 
missing momentum p m [MeV/c]. The solid line corresponds to the RH, the dash-dotted line to the 
RHF and the dashed line to the RDDH approach. 
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FIG. 2. Left panel: density distribution r 2 n p i/2(r) in [/m _1 ] as a function of the radius r ([fm]). 
Right panel: density distribution in momentum space n pl / 2 (Pm) in [GeV/c~ 3 ] as a function of the 
missing momentum. The lines are denoted as in Fig. 
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FIG. 3. Left Panel: upper component g p \j2 (Pm) of the relativistic bound state wave function in 
[(GeV/c)~ 3 / 2 ] as a function of the missing momentum. Right panel: nonrelativistic limit Qpl^ip)- 
The lines are denoted as in Fig. |l[ 




FIG. 4. Left panel: lower component f p i/2(Pm) of the relativistic bound state wave function in 
[(GeV/c)~ 3 / 2 ] as a function of the missing momentum. Right panel: nonrelativistic limit fpl^ip)- 
The lines are denoted as in Fig. |l[ 
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